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1. Historical 


Conformal mapping has been used as a tool for solving 
problems in fluid mechanics and electromagnetism for more than 
one hundred years. Riemann's (somewhat incomplete) proof of the 
possibility of mapping closed contours on circles dates from 1850 
[1]. Schwarz introduced what is now known as the Schwarz- 
Christoffel transformation in 1869 [23. In the same year, Kir- 
chhoff and Helmholtz used conformal mapping to solve classical 
problems of flows with free surfaces [33. The importance of con- 
formal mapping in fluid mechanics in the second half of the 
nineteenth century and the first quarter of the present one stems 
from the well-known property of invariance of the Laplace equa- 
tion through conformal mapping; the mapping itself is, indeed, 
defined by a Laplace equation. Since the theory of images allows 
any incompressible, potential flow defined by given singularities 
to be determined in the presence of a circle or a straight line, 
the analogous problem in the presence of an arbitrarily given 
contour is reduced to the problem of mapping the contour onto a 
circle or a straight line. In so doing, the problem of determin- 
ing the analytic function representing the complex flow potential 
is split into two parts. One consists of finding the potential 
for the circle, and that is easy. The other still requires the 
determination of an analytic function to define the mapping, 
which is as difficult a task as the original problem [4, page 
773. Nevertheless, it is much more appealing to people who have 
a visually oriented mind, and it can be subdivided in a number of 
successive steps, each one of which is easy to understand 
geometrically. 


The Joukowski mapping, 


discovered in 1910 [53, is extremely simple and easy to handle. 
Using a circle as the basic contour in the ^-plane, it can gen- 
erate a two-parameter family of cusped profiles and a three- 
parameter family of smooth contours ranging from familiar el- 
lipses to peanut-shaped figures (Fig. 1). Despite their poor 
aerodynamical properties, the Joukowski profiles played a crucial 
role in the understanding of the mechanism of lift and, conse- 
quently, in the establishment of the theory of flight as a sci- 
ence. In creating his profiles, Joukowski introduced an idea 
which turned out to be extremely rich in consequences: that is, 
that a shifting of the center of the circle is sufficient to pro- 
duce an airfoil having either camber or thickness or a combina- 
tion of both. 
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form: 


Another major contribution dates from 1918 when von 
and Trefftz [6] observed that, if (1) is recast in the 
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it can be interpreted as a combination of two bilinear mappings 
and a power. If 6 is used instead of 2 in (2), the KArmSn- 
Trefftz mapping results: 
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which, for 6=2-e/ir, produces an airfoil with a finite internal 
angle, e, at the trailing edge. The property of a power, 

z = c 6 (4) 


to eliminate sharp corners, was known and had been applied since 
the times of Schwarz and Christoffel, but the importance of its 
application in (3) should not be underestimated. 


The role played by the exponential (and, conversely, the 
logarithm) and by all functions directly related to it (direct 
and inverse trigonometric and hyperbolic functions) in generating 
infinite cascades can be traced back to Kflnig, 1922 [73. 

Transformations of this kind, for example, 
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have essential singularities at infinity. Therefore, the image 
of the physical point at infinity in the plane of the circle is 
generally a pair of spiral-vortices at a finite distance, a 
geometrical property related to the ability of a cascade to de- 
flect the flow. 

Finally, in 1933 Ferrari found the basic mapping for bi- 
planes. He showed that a doubly periodic function was needed, 
that is, an elliptic function. His paper [8] is the first appli- 
cation of elliptic functions to wing theory. [8]. 

Important as these basic mappings were, they did not 
solve the problem of mapping an arbitrarily given contour onto a 
circle. They generate a closed contour without sharp corners 
which, in a favorable case, may look like a quasi-circular con- 
tour which reminds one of a potato (the word "potato" has actual- 
ly been used occasionally in the literature to designate such 
figures and, although not defined in mathematical terms, seems to 
be as acceptable as "quasi-circle"). The technique for the map- 
ping of the quasi-circle onto a circle, proposed by Theodorsen in 
1931 [93, actually brought conformal mapping from the speculative 
level down to practical levels of aeronautical engineering. Its 
basic idea is indeed extremely simple; if the center of the cir- 
cle and the centroid of the quasi-circle are located in the ori- 
gins of their respective planes, and some scaling is used to make 
the areas of the two figures coincide, then the mapping can be 
expressed in the form: 


X 4) i 0 

where z=re \ c=e and the modulus of f is small. Therefore, 
two equations follow: 

In r = A ( 0 ) f $-0 = B(e) (7) 

if f = A + iB. On the circle, A and B are conjugate Fourier 

series. If In r were known as a function of 6, A could be deter- 

mined by expanding In r into a Fourier series, and B would follow 
automatically. Since In r is known, instead, as a function of 4> , 
iterations are necessary to make the two equations compatible. 

Alternatively, values of $-e can be obtained rephrasing the prob- 

lem as an integral equation, so that the formula to be solved by 
iteration is the Poisson integral: 
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dw 


Convergence of the procedure is assured if the potato is "star- 
shaped”, that is, if its contour is crossed only once by any ra- 
dius issued from the center [10]; therefore, the procedure is 
generally safe and efficient. 

Applications of the Theodorsen method to airfoils, bi- 
planes and cascades were pioneered by Garrick [11 t 12,131* who 
used the transformations introduced by Joukowski, Ferrari and 
Kflnig to produce quasi-circles as intermediate images of given 
contour s. 

For a better understanding of the rest of this paper, 
three points should be made. 

1) Happing of a contour into an _exa_c_t circle was neces- 
sary, when problems of incompressible, potential flow had to be 
solved, in order to take full advantage of the invariance of the 
Laplace equation and of the easy evaluation of the flow past a 
circle, 

2) Since the only values of interest belonged to the ri- 
gid contours, all calculations could be limited to the circumfer- 
ence of the circle. 

3) Nevertheless, all calculations, with the exception of 
the ones related to the Joukowski mapping, were extremely cumber- 
some when performed by hand. Host of the time, the Joukowski 
mapping was preferred to the Karman-Tref f tz mapping. Cascade and 
biplane problems were rather analyzed by approximate methods. 
The Hilbert integral, with the entire circumference subdivided in 
just a dozen intervals, was preferred to the more laborious 
Fourier expansions. 

Such practical difficulties and the shifting of interest 
towards problems in compressible flow contributed to relegate 
conformal mapping into the background in the forties. When com- 
puters came about, there was not so much a demand for solution of 
incompressible, potential flow problems, and the consequent ap- 
plication of conformal mapping techniques. 
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2 . Conformal mapping as a grid generator 


Nevertheless, in recent times conformal mapping has again 
been looked upon in connection with problems of compressible 
flow, but of course the use of it as a device to solve the same 
equations of motion in a simpler plane has been dropped, since 
neither the Euler equations nor the Navier-Stokes equations are 
invariant by conformal mapping. The reason for revival stems 
from the fact that in contemporary numerical gas dynamics the 
equations of motion must be solved iji the entire flow field , be- 
cause of compressibility effects, viscous effects, presence of 
vorticity, etc. Contrariwise, for incompressible, potential flow 
it is theoretically correct and practically convenient to search 
for the solution on the rigid contour only . A computational grid 
is necessary, upon which the equations are to be discretized. 
When two space-like variables are involved, conformal mapping is 
a very convenient tool to generate a computational grid. 

If the boundaries of the flow field are mapped onto a 

i A 

circle centered at the origin of a complex plane s = pe , the 
network of p=constant and 6=constant lines is conversely mapped 
onto the physical plane as an orthogonal grid t well draped around 
the rigid contours because the latter are the image of the circle 
itself. 


Since the grid is orthogonal, the vector operators in the 
physical plane can be expressed directly and easily in terms of 
their counterparts in the mapped plane. 

The search for an appropriate mapping is actually a 
search for a single function of a complex variable, a task by far 
simpler than the search for two functions of two independent 
variables. All the findings of more than a century, properly 
digested and interpreted as briefly shown in Section 1, can be 
put together to suggest the most appropriate mapping for a par- 
ticular task. Our computers have eliminated the tediousness, 
inaccuracy and material errors which used to hamper the hand cal- 
culations of half a century ago; they also allow the necessary 


data for a very refined grid to be evaluated in a fraction of a 
second. Complex arithmetic in FORTRAN and other languages of the 
same level reduces coding to just about a rewriting of the basic 
equations (except when angles contained in more than one quadrant 
are involved). 

Practical applications of conformal mappings to generate 
grids have been developed in the last decade following two 
separate lines. One proceeds from the observation that any con- 
formal mapping is defined by an analytic function, and that real 
and imaginary parts of the latter are harmonic functions. There- 
fore, the conformal mapping problem consists of solving two La- 
place equations. The task can be formulated in strictly numeri- 
cal terms, using modern, fast Laplace solvers. This viewpoint 
can be traced back to 1923 [14] and has been made popular 
by a 1975 paper by Thompson et al. [15] . The other puts the em- 
phasis on the use of closed form analytical expressions for the 
mapping functions. It seems that the first non-airfoil related 
application of the technique was presented in 1972 [16], although 
the same problem was reformulated in a more general form only in 
1974 [17]. We will try to analyze here the philosophy of the ap- 
proach and to show how it works in practical cases. 


3 . Philosophy of the ' closed form 1 approach 


We will begin by discussing some of the advantages of a 
closed form approach. To fully appreciate them, let us focus our 
attention on what the majority of problems of current interest 
are: Two-dimensional, unsteady flows or three-dimensional, steady 
flows. In the former, the physical space is two-dimensional but 
the computational grid may be variable in time; in these cases we 
need a grid, generated by conformal mapping at every instant of 
time, but depending on parameters which are functions of time. 
In the latter, we may find it convenient to create a computational 
grid on planes defined by two of the three space coordinates, and 
again letting the grid change as a function of the third coordi- 
nate. To organize the following formulae, let us stipulate that 
the Cartesian coordinates on any plane to be mapped are x and y, 
and that a complex coordinate is defined, 



z s x + i y 


( 9 ) 


using the symbol t either for time or for the third space coordi- 
nate. The mapping of each z-plane onto another complex plane, 
defined by a variable £, will then be accomplished by an analytic 
function : 

5 = c(z;c> (10) 
where c is shorthand for any parameter, function of t. 

As in every flow problem, it is not so much the coordi- 
nate of a point as the derivatives at that point which matter. 
First, we need the derivatives of the coordinates in the physical 
plane with respect to the coordinates in the mapped plane, and 
vice versa, contained in the complex expression: 


g 



( 11 ) 


Then, we need second derivatives, which appear in dealing with 
the curl of the velocity vector, and these entail dg/dz. Final- 
ly, we need derivatives with respect to t, viz. 3e/3t and 3g/3t, 
which are computed from (10) and (11) by differentiating with 
respect to c and keeping z constant. Obviously, it is very con- 
venient to have the mapping expressed in a closed analytical form 
since all derivatives are also expressed in closed form and can 
be exactly evaluated where and only where they are needed, that 
is, at the computational nodes. The argument is particularly in- 
teresting with regards to the t-derivatives. If the grid changes 
in t, the values of z generally change in t at each nodal point; 
therefore, numerical differentiation at a constant z may be 
cumbersome and inaccurate. 


There is another case where having a closed form mapping 
is convenient. Potential (but not incompressible) flow problems, 
with the flow field extending to infinity, are commonly solved by 
relaxation. Commonly, the physical potential is expressed as the 
sum of an unknown and the potential of a flow, satisfying the 
proper physical conditions at infinity. The latter is easily 
formulated in terms of the transformed variables if the mapping 
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is defined by a closed form function. 


Then, a difficulty seems to arise. As we said in Section 
1, in most cases it is not possible to find a closed form mapping 
which provides an exact circle as the image of a given contour; 
it is clear that Theodorsen's step is not a closed form function. 
Nevertheless, the difficulty is only apparent. Let us assume, 
indeed, that we know a function capable of transforming a given 
contour into a quasi-circular potato. Let us consider a computa- 
tional plane with two variables, X and Y, relating X to p and Y 
to 9. The functions X(p) and Y ( 0 ) will be so defined that X=0, 
Y=0, X=1 and Y=1 on each boundary, in turn. The computational 
region in the (X, Y) plane is, thus, the interior of a square. 
In most cases, the functions X(p) and Y(6) are used to apply 
proper stretching of coordinates in either direction, in order to 
concentrate nodes where necessary. The computational grid is 
orthogonal and divided into equal intervals. The corresponding 
grid in the c -plane is not orthogonal (not only because the X=0 
line is not an exact circle, but also because the other boun- 
daries may not be circles or straight lines either). Consequent- 
ly, the grid in the physical plane is not orthogonal. Accuracy, 
however, is not impaired for want of orthogonality. The equa- 
tions, originally written in the (p, 0) frame, are recast using X 
and Y as independent variables. Some additional terms will ap- 
pear; the very important boundary condition on the rigid surface 
will be properly written by stating that the normal velocity com- 
ponent (not just the p-component) vanishes. In conclusion, we 
are not expecting any dramatic advantage from having a circle as 
one of the boundaries, and therefore we may consider our mapping 
problem solved when we find a quasi-circle as the mapped image of 
the given boundary. The Theodorsen step can be dropped with all 
its additional burden of iterations, Fourier expansions, spline 
fittings , etc . 

Let it be clearly stated that we are not afraid of unsur- 
mountable difficulties or unaffordable computational times con- 
nected with the Theodorsen step. That old (but hard to execute 
by hand) way of computing Fourier expansions, which has been 
given a new popularity under the FFT label serves the purpose 
egregiously well. Dropping the Theodorsen step, whenever possi- 
ble, is justified by our desire of achieving a solution to the 
problem in a closed form, with a view to the formal calculation 


of derivatives, particularly with respect to t. There is a de- 
finitive need for them in all problems where the grid depends on 

t . 


This philosophy was clearly exposed in [16] and accepted 
by Jameson [18,191. We will discuss some applications in Sec- 
tions 5 through 8, Before that, we will mention some of the 
techniques which we tend to classify as numerical, rather than 
analytical . 


Conformal mapping as a Neumann problem 


The closed form approach may look like an empirical at- 
tempt to solve the problem, on the basis of analogies and ima- 
ginative variations, and therefore strongly dependent on the men- 
tal structure of the investigator. I have heard the word "art" 
used in a derogatory sense in connection with this type of work, 
on other occasions. Actually, good science is the product of 
personal ingenuity and crafty skill. What tends to be classified 
as "scientific" in these days is rather "technological", that is, 
some process which has been sealed in a black box for general 
purposes . 


Whatever the consequences, there is no doubt that the 
quest for an organized conformal mapping procedure is legitimate. 
It is also classical. Green’s formula: 

= ° < 12 > 

for two harmonic functions, u and v, regular in the domain sur- 
rounded by the closed contour on which the integral is made, and 

on the contour itself, dates from 1828. With any two points on 

i 0 i<b i a 

the contour denoted by £=pe and z=re , and with z-z =ae , it 
follows from (12) that 

u(z) = - 1 j [uU) in a] ds ( 13 ) 


All formulae related to conformal mappings can be obtained from 
(13); 3u/3n is generally known from its conjugate, the variation 


of the tangent along the contour. Then (13) becomes an integral 
equation for u. 

Many different forms of the equation can be obtained, if 
one makes use of well-known properties and formal rules, such as: 
the Cauchy-Riemann conditions, integration by parts, Schwarz and 
Poisson's integrals, integrals defining the coefficients of a 
Laurent series. Different forms are also obtained by taking the 
basic contour as a circle or as a straight line, and by defining 
u either as the logarithm of z or as the logarithm of g. For ex- 
ample, if we start from (6) and from a circle, the other contour 
being a quasi-circle, we obtain the Theodorsen mapping in its in- 
tegral form (8); the Fourier series form follows easily. The same 
equation (6) and somewhat different integral equations have been 
used by Symm [20] and Hayes et al . [21] to produce numerical 
techniques which are not restricted to mappings of quasi-circles 
onto circles, as in Theodorsen's, but apparently can handle any 
(probably, star-shaped) contour. If we use the logarithm of g, 
we can interpret Theodorsen's ideas in terms of derivatives, a la 
Timman [22], a variation which seems to offer some numerical ad- 
vantages [17]. If we start from a straight line and an arbi- 
trary, closed contour, using again the logarithm of g, we obtain 
what Davis [23] presents as a generalization of the Schwarz- 
Christoffel mapping to a polygon with an infinite number of in- 
finitesimal vertices: 

log g = | / log (c-b) dB (1*0 

Anyone interested in these types of comparisons could profitably 
read a paper by Birkhoff et al . [2*1] which is outdated only from 
a computational viewpoint. 

Proceeding in the opposite direction as Davis, the 
Schwarz-Christoffel formula for a polygon with a finite number of 
vertices can be found as a particular case of (1*1). This formula 
is, in principle, a very powerful mapping tool. It can map a po- 
lygon on a circle, without restricting number, location and aper- 
ture of vertices, or the lengths of the sides, and permitting 
vertices to be located at infinity. It is really a definition of 
g, rather than of ?, which must be obtained by complex integra- 
tion in a numerical form, in almost all the cases. This is not a 
major shortcoming, however, since numerical integration can be 


performed quickly and efficiently. The derivative, dg/dz is 
straightforward. Nevertheless, it is known that the coeffi- 
cients, ^ and which appear in the formula 


g = n U - 5 .) 1 (12) 

must be obtained by trial-and-error iterative processes (see, for 
example, [23,25,26]). If the grid does not depend on t, its 
coefficients can be determined once and for all. In this case 
the Schwarz-Christoffel formula belongs to the category which we 
consider in the present paper. It does not if the grid depends 
on t. The same may be said for all mappings obtained by solving 
a Neumann problem via iterations on an integral equation or 
Fourier expansions. 


5 . Karman-Tref ftz mapping for airplane cross sections 


We will now consider some mappings using a finite number 
of closed form relations. 

The numerical analysis of a steady, supersonic flow past 
an airplane may be performed by marching in an axial direction 
and updating values at successive cross-sectional planes. The 
region of interest in each plane is bounded by the section of the 
body and the section of the bow shock. Conformal mapping of the 
body onto a quasi-circular shape provides a grid which tends to 
become a polar grid at infinity and therefore is the best suited 
to adjust to the shape of the body, whatever it is, and to the 
almost circular shape of the bow shock. The body shows a number 
of edges and corners as those indicated by letters in Fig. 2. If 
the airplane is arrow-winged, stations will be reached where the 
body will be composed of three unconnected parts; if two ficti- 
tious lines are drawn between the fuselage and the trailing edges 
of the wings, again we can see corners and edges at all points 
denoted by letters in Fig. 3. Observing that the corners and 
edges always come in pairs, because of the symmetry of the cross- 
section, we can think of eliminating them by successive applica- 
tions of the Karman-Trefftz mapping (3), with the singular points 


either on the contour or, in case of rounded edges and corners, 
slightly inside. Note that edges require values of 6 between 1 
and 2 and corners require values of 6 less than 1. When the 
latter are applied, the corner is open, but so is the rest of the 
plane, part of which may end up in a second Riemann sheet. In 
principle, this is not an obstacle to the removal of corners, be- 
cause the portion of plane which disappeared will be recalled 
when removing the next edge. In practice*, a quite cumbersome ad- 
ditional piece of logic must be added to identify points belong- 
ing to the second Riemann sheet. The trouble can be avoided by 
executing the mappings not according to the order of appearance 
of a corner or edge along the contour but in a sequence of de- 
creasing values of 6. 

From the viewpoint of coding, the repeated application of 
the Karman-Tref f tz procedure has many advantages: 

1) Regardless of the number, position and aperture of corners 
and edges, the same operation is used, which means the code may 
be written in the form of a loop and applied as many times as 
necessary , automatically, 

2) The mapping can be inverted, and the inverse mapping has the 
same form; therefore, the same routine can be used for the direct 
and the inverse mapping, 

3) The derivatives are easily coded; for example g, as defined 
by (11), is actually the product of the derivatives of each in- 
termediate step, and the logarithmic derivative of g is the sum 
of the logarithmic derivatives of such steps. 

To show how close to a circle the image of a fuselage 
with two sections of arrow wings is. Fig. presents a computa- 
tional grid in the physical plane and its image in the mapped 
plane. Naturally, with the bow shock very close to the leading 
edge of the wings, its own image is far from a circular shape 
but, as we said above, departure from orthogonality of the grid 
is not jeopardizing accuracy, and this is particularly true in 
the vicinity of the bow shock, where the flow is uncomplicated. 
Details of the technique and its application to the arrow-winged 
airplane problem can be found in [17]. 


6. Imaginative devices 


There are no limits to the number of shapes which can be 
obtained by executing elementary mappings in a sequence and using 
a little ingenuity. Indeed, bilinear transformations, powers and 
logarithms are the building blocks with which one should learn to 
play, always keeping in mind that rotations and translations can 
be cleverly used to locate singularities where needed. Here is 
an interesting example, due to Rossow [27], which is entirely ex- 
pressible by a sequence of bilinear mappings and powers (Fig. 5). 

A circle, centered at the origin, is translated upwards 
and then a Joukowski mapping changes it into an arc of a circle, 
counted twice. The arc is rotated about one of its ends and a 
new singular point is defined somewhere along its length. A new 
Joukowski mapping is applied, in reverse, so that the portion of 
the arc between the two singular points become a circle again, 
and the portion left outside remains appended like an infinitely 
thin tail. Finally the circle is relocated, and a third 
Joukowski mapping is used to transform it into a Joukowski pro- 
file; the little tail becomes a flap or spoiler, whose location 
and length can be controlled by changing the parameters used in 
the successive steps. According to our quasi-circular philoso- 
phy, this very simple mapping can be used for any airfoil with 
attached (but not infinitely thin) spoilers or flaps. 

I faced a similar problem when confronted with generating 
a grid for the calculation of the precursor muzzle blast [28]. I 
needed a grid shaped as in Fig. 6; the contour defined by D'C'CD 
can be mapped onto the real axis of a w-plane, and its exterior 
onto the upper half w-plane by a simple Schwarz-Christoffel 
transformation [29, page 159], but this mapping would not provide 
a family of grid lines issuing from what has to be interpreted as 
the bore of the gun and wrap around the barrel, as in the figure. 
The problem was solved by defining two new singular points, B and 
B', in the w-plane and applying to this plane an inverse 
Joukowski mapping onto a £-plane. Radial straight lines and con- 
centric circles in the s-plane are now producing the wanted grid 


in the z-plane. Note that the position of B and B', being arbi- 
trary, permits the ratio of outer-to-inner radius of the barrel 
to be matched (Fig. 7). With the change of notation: w = Z + 
1/Z, which simplifies the coding, the mapping is thus defined by: 


z = (r /ir) [ (Z 2 -1/Z 2 )/2 = log 2 ? - iirj 
Z + 1/Z = 2 B (c + I/O 


(15) 


In the problem just described, the computational region 
is limited, on one side, by the precursor shock which moves out 
in time. Therefore, the computational grid is a function of 
time, but the dependence on time shows only through the stretch- 
ing parameters; the mapping itself remains invariable in time and 
the equations of motion carry no terms of the type 3g/3t or 
3c/3t. When a mapping is needed for the same problem in the 
presence of a protruding bullet, however, the contour itself 
changes in time and so does the mapping. Once more, the problem 
can be solved with little additional effort. One can start with 
a half circle (Fig. 8), reduce it to a half ellipse using a 
Joukowski mapping, and then apply a Karman-Trefftz formula to 
change the angles between the contour and the real axis. After 
that, the mapping continues as defined by (15). The axis ratio 
of the ellipse keeps growing as the bullet nose advances. The 
power in the Karman-Trefftz mapping decreases from an initial 1 
to a minimum value of 1/2. At this stage, the contour in the 
physical plane makes a 90 degree turn at B and B', as required to 
accomodate the grid to the side surface of the bullet (Fig. 9). 


]_, Different mappings for the same geometry 


In solving problems of incompressible, potential flow 
there is no ambiguity about the choice of the mapping since the 
contours of interest are specified exactly and the mapping func- 
tion has to be regular over the entire flow field (boundaries ex- 
cluded). In the present context, though, different mappings may 
accomodate the same contours but generate completely different 
grids, as we have mentioned in the preceding Section. Care must 
be taken to use a mapping whose grid is the best suited for the 


problem in hand. This problem is currently acute in cascade 
analysis, as we will mention later; here I would like to show it 
in a simpler connection. 

We want to generate a grid for the intake of which Fig. 
10 shows the centerline and the shroud. We can drape our grid 
about a semi-infinite slit, parallel to the real axis in the 
upper half-plane. Two mapping functions, apparently opposed to 
one another, can serve the purpose. The first is: 

z = 5 + e s (16) 

This function maps the slotted half plane onto a strip, which is 
a very convenient domain for our computational variables X and Y. 
The corresponding grid (Fig. 11 ) is convenient for the interior 
of the intake, but it needs some stretching to provide resolution 
to the exterior of it; in particular, the outer surface of the 
shroud is poorly resolved. The other mapping is defined by: 

z =e + a log 5 ( 17 ) 

and it maps the same region of the physical plane onto the entire 
upper c-half-plane . Cartesian coordinates in the s-plane produce 
the grid of Fig. 12. In this case, the grid is very good outside 
and very poor inside. If we start from polar coordinates in the 
C-plane, we obtain the grid of Fig. 13. The general appearance 
of the grid lines recalls Fig. 11, but the situation is reversed: 
the resolution is very poor inside and very good outside; in this 
case, an accumulation of p=constant lines near the origin of the 
S-plane is necessary to generate some coordinate lines inside 
the intake. An application of (16), with a stretching of coordi- 
nates to relax the resolution inside the nacelle, has been made 
by Caughey and Jameson [30] . 

Anyway, these mappings are obtained in a very straight- 
forward manner, and more complicated manipulations [31] do not 
seem necessary. 


£3. Biplanes, revisited 


Ives [32] has proposed a technique for the mapping of two 
airfoils which is rich in possible consequences. We have seen in 
Section 5 how the successive application of Karman-Trefftz map- 
pings can eliminate corners and edges on a contour, producing a 
quasi-circular shape. A similar idea, exploiting the theory of 
images with respect to a circle, allows two airfoils to be mapped 
on two concentric circles. Two Theodorsen mappings are used as 
intermediate steps, as follows (Fig. 1*0. First, a Karman- 
Trefftz formula is applied to transform one of the airfoils into 
a quasi-circle, and then the quasi-circle is transformed into an 
exact circle by the Theodorsen technique. At this stage, the 
second airfoil is still shaped as an airfoil although with a dif- 
ferent shape. The next step manages to transform the second air- 
foil into a quasi-circle, without distorting the first circle . 
The problem is solved by using the product, side-by-side, of two 
Karman-Trefftz formulae; the first contains the two singular 
points pertinent to the second airfoil (a and b) and their coun- 
terparts in the mapped plane (a and $); the second contains the 

## 2 # 2 # 

images of such points (1/a , 1/b and r /a ,r /0 , respectively, 

O V 

where conjugates are denoted by * and r Q is the radius the image 
of the first circle in the mapped plane): 


(C-a)U-r 2 /a*) 

(C-0)U-r 2 /0*) 

o 


j- (z-a) (z-1/a ) j<5 
(z-b) (z-1/b*) 


(18) 


Finally, the second quasi-circle is moved inside the first circle 
by a bilinear transformation in such a way that its centroid 
coincides with the center of the circle, and a second mapping of 
the Theodorsen type is applied. Again, such a mapping must take 
into account that the flow region is a ring between two concen- 
tric circles. Therefore, the exponent in (6) cannot be a simple 
Taylor series, as it would be if the flow occurred inside a cir- 
cle, or a Taylor series of negative powers, as it would be if the 
flow occurred outside a circle, but it must be a two-sided 
Laurent series; we express this need by saying that the Theodor- 


sen mapping is defined in this case by: 


z = ? g f (c ;c)-f( i/t ;c ) ( 19 ) 

where c is any constant occurring in f(?) [ 12 ], 

The idea of using the images of the singularities with 
respect to a circle in order to operate on the other contour 
without distorting the circle can surely be extended to problems 
dealing with two separate contours but not necessarily airfoils. 

Two questions can be posed. First, how can Ives' tech- 
nique produce the mapping, without using elliptic functions? The 
explanation can be found in the fact that the sequence of map- 
pings used by Ferrari and Garrick consists of a logarithm, fol- 
lowed by an elliptic function. The net result actually has only 
one period, the other being neutralized by the multi-valuedness 
of the logarithm. On the other hand, the periodicity is intro- 
duced in Ives' mapping through the play of reflections of singu- 
larities produced by (19). The second question is whether one 
could bypass the two Theodorsen corrections, in the spirit of 
Section 3 above. Of course, in this case we would have to deal 
not with one distorted circle alone, but two, the first being 
particularly important since the images of ( 19 ) would, in any 
case, be defined with respect to a circle which now would only be 
an approximation to the real contour. I believe that the tech- 
nique could still be applied, but no examples are available. 


9. Cascades 


The climax of difficulties is reached in the problem of 
generating a grid to compute flows through a cascade of airfoils. 
Garrick [131 sensed the difficulty well ahead of the computer 
era, but he did not have a way of measuring it; as we said, hand 
computations were necessarily limited to the simplest cases. He 
said: "It is to be noticed that improvements in the initial 
transformation are desirable and should be sought, particularly 
to take care of highly cambered airfoils more conveniently, in 
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order to reduce the amount of subsequent calculations." At that 
time, the initial transformation was the one defined by (5), ob- 
viously a poor way to get started when dealing when highly cam- 
bered airfoils, since the basic shape furnished by (5) is a cas- 
cade of flat , straight, double-sided segments. Unfortunately, in 
a cascade four parameters must be considered, solidity, stagger, 
camber and thickness, and these parameters interact with each 
other, whatever the choice of the basic mapping. So long as 
solidity and stagger are low, even ( 5 ), applied to a moderately 
thick and cambered airfoil, produces a reasonable quasi-circle. 
But as solidity and stagger increase, the contour tends to become 
peanut-shaped (Fig. 15), with a catastrophic distribution of 
points around it (something like a wide circle, corresponding to 
a small portion of the original profile, with a small appendix, 
which is the image of all the rest); and, at times, the contour 
is not even star-shaped any longer. Can (5) be blamed for such a 
behavior? A detailed discussion of this question would transcend 
the limits of the present paper; let it just be said that the 
difficulty does not disappear when another mapping function, 
adopted by Legendre [33] and Ives [34], is used (Ives' function 
is the same as Legendre, to within a rotation and a bilinear 
transformation). This mapping can be written in the form; 

r ? + b-,5 sin A(z - 1 ) fofn 

L C - b J " sin A( z + 1) J 

where z and c are the physical plane and the mapped plane, 
respectively, 1 and -1 are the location of the trailing edge and 
the center of the leading edge of a profile, respectively, 6 is 
the usual Karman-Trefftz exponent, A is related to solidity and 
stagger by 


2 i A = tt s e lB (21) 

and 

b = - i tan 7 ( 22 ) 

So far, the major advantages of this transformation with respect 
to ( 5 ) are the presence of an exponent, 6 to produce profiles 
with a finite trailing edge angle, and the possibility of comput- 
ing either c as a function of 2, directly from (20), or 2 as a 
function of £ from the inverse of (20): 


z 


( 23 ) 


, 1 , e 2i# -w 8 

= - 1 - 2A 106 ~^U s 

e - w 


with 


w 


S ± b 
c - b 


(2U) 


When stagger and solidity are high, most of the skill which must 
be used to follow the contour properly in applying (20) is wast- 
ed, due to the ungainly shape of the resulting contour. More 
manipulations are needed before a quasi-circular shape is ob- 
tained; but we will not elaborate on this point in the present 
paper. 

Let us assume then, that a circle has been obtained 
somehow from the given cascade; and, in the spirit of Section 3, 
we will gladly accept a quasi-circle, without refining the map- 
ping any further. At this stage, we must decide what type of 
grid we want. Different choices can indeed be made for the gen- 
eral appearance of the grid. Three of these are the most common 
(Fig. 16). The first contains lines which run from left to right 
through the cascade, more or less as streamlines of a real flow 
would do. The second (commonly called an 0-grid) contains lines 
wrapped around each profile. The third (commonly called a C- 
grid) contains similar lines which run around the profile and an 
infinite line issuing from the trailing edge, as a wake. Now, 
regardless of the choice of the first mapping step, the plane of 
the circle will always have four singular points, representing 
the points at infinity (on the left and the right of the cascade) 
and their reflections on the circle. If we invert the circle 
plane so that the interior of the circle corresponds to the 
exterior of the cascade, then the two singular points inside the 
circle represent the physical points at infinity. A system of 
streamlines and equipotential lines for an incompressible, poten- 
tial flow, proceeding from one singular point to the other and 
having stagnation points at the points on the circumference which 
correspond to the leading and trailing edge, provides a grid of 
the first type. In building it up, one finds numerical difficul- 
ties. It is indeed hard to follow a streamline on a spiral vor- 
tex when the singular point is very close to the periphery of the 
circle (this is what happens with high solidity and high 
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stagger) . 


Polar coordinates are very simple to use, but they gen- 
erate a very inconvenient grid because of the presence of a 
singularity in the flow field between profiles, at the point 
corresponding to the center of the circle [35]. (See Fig. 17.) 

For an O-grid, the best procedure has been suggested by 
Ives [35]. It consists of mapping the circle onto a rectangle, 
in such a way that two vertices correspond to the singular points 
(points at infinity in the physical plane) and the opposite side 
corresponds to the circumference (the contour of the profile). 
The basic function defines the correspondence shown in Fig. 18 
(where there is a straight slit between the two singularities). 
In the same figure, lines corresponding to straight lines, paral- 
lel to the sides of the rectangle, are shown inside the circle. 
The function which performs the task is the simplest of the Jaco- 
bian elliptic functions: 


z = sn (e , m) (25) 

where m is defined by the position of the singularities and tends 
to 1 when solidity and stagger increase. Using Landen's 
transformations, the sine-amplitude can be expressed in terms of 
trigonometric functions (for m close to 0) or the hyperbolic 
tangent (for m close to 1); the coding of the subroutine is obvi- 
ous and the computational time is negligible. Therefore, (25) is 
a convenient function. 

Note that the circle (or quasi-circle) obtained through 
(5) or (20), and the circle obtained from (25) are by no means 
the same, despite the fact that in both planes the images of the 
points at infinity are symmetrically located on the real axis. 
The second circle is centered at the origin, but the first is 
not; therefore, a bilinear transformation must be used to map the 
two circles onto each other. This mapping depends on all four 
basic parameters and the trailing edge angle as well. Its in- 
terpretation, thus, is not easy, but a systematic study is needed 
to understand what range of basic contours, acceptable as a back- 
ground of a grid for given profiles, can be obtained by using 
(25), a bilinear mapping and (23), in that sequence, without 
resorting to a Theodorsen step. It seems to me that we either 


have reasonable shapes, in which case that step can be skipped, 
or combinations of solidity and stagger which tend to generate 
intermediate shapes similar to the one on the right hand side of 
Fig. 15, in which case the Theodorsen technique would not work, 
and where another intermediate step, of a different nature, 
should be used. Even Joukowski mappings may help (see the last 
contour in Fig. 1). Work along these lines is urgently needed, 
if we want to obtain simple mapping procedures for three- 
dimensional turbomachinery problems. To make the point, we 
present a three-step sequence, in which the first circle, ob- 
tained from the rectangle, is the one of Fig. 18. A second cir- 
cle (Fig. 19) is obtained by a bilinear transformation and the 
cascade, obtained through (23) is shown in Fig. 20. The contour 
is ugly, but the grid is perfectly usable. 


10. Grids for ablated , three-dimensional bodies 


We conclude this presentation showing a method applied to 
generate a computational grid for a three-dimensional, time- 
dependent problem. The flow to be determined is the shock layer 
around the ablated nose of a cone-cylinder; the flow is mostly 
supersonic, but it may have a subsonic bubble and an imbedded 
shock. The geometry of the body, which is ax isymmetrical before 
ablation, becomes three-dimensional because of different ablation 
in different meridional planes. A grid is needed in a number of 
these planes, and corresponding points must be connected between 
adjacent planes, in order to generate a three-dimensional compu- 
tational mesh. It is therefore convenient to have a grid defined 
in closed form, to make the evaluation of circumferential deriva- 
tives as easy as possible. For a given section of the body (Fig. 
21) we define a skeleton, that is, a polygon, all contained in- 
side the body, approximating the shape of the wall. Instead of 
using a Schwarz-Christof f el function to transform the skeleton 
into a straight line in a single operation, we have opted for 
straightening one vertex of the polygon at a time, beginning with 
the one farther from the nosetip. As in the case of the repeated 
Karman-Tref ftz mappings, the procedure is easily coded in a loop. 
The last step consists of a square root transformation, in order 
to bring the skeleton on the real axis and the body axis on the 


imaginary axis of the ?-plane. The image of the wall on the 
C-plane is close to a straight line. The grid, which will be 
normalized between the image of the wall and the image of the bow 
shock, is always well shaped, despite strong concavities of the 
body produced by severe ablation. The method, first tested on 
axisymmetric problems [37], has been successfully applied to 
three-dimensional problems [ 38 ]. 
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Figure 7,- Generation of the muzzle mapping 



Figure 8.- Generation of the muzzle mapping 
with protruding bullet. 



Figure 9.- Grid for muzzle with 
protruding bullet. 



Figure 10.- Nacelle geometry. 



Figure 


11.- Grid generated by equation (16). 



Figure 12.- Grid generated by equation (17), 
Cartesian coordinates. 






Figure 17.- Frith's grid for cascades. 



Figure 18.- Circle mapped from rectangle, 
equation (25) . 
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Figure 19.- Circle mapped from cascade, equation (20). 



Figure 20.- Cascade obtained from circle of Figure 19. 



